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Abstract 

We present a first scaling test of twisted mass QCD with pure Wil- 
son quarks for a twisting angle of tt/2. We have computed the vector 
meson mass my and the pseudoscalar decay constant Fps for differ- 
ent values of (5 at fixed value of romps- The results obtained in the 
quenched approximation are compared with data for pure Wilson and 
non-perturbatively 0(a) improved Wilson computations. We show 
that our results from Wilson twisted mass QCD show clearly reduced 
lattice spacing errors, consistent with 0(a) improvement and without 
the need of any improvement terms added. These results thus provide 
numerical evidence of the prediction in ref. [1]. 



1 



1 Introduction 



Formulating QCD on a space-time lattice admits a substantial amount of free- 
dom in discretising the continuous derivative with restrictions coming from 
obeying principles such as gauge invariance, locality and unitarity. Naturally, 
as long as universality holds, all these formulations should provide consistent 
results in the limit when the discretisation is removed. The standard Wilson 
formulation of lattice QCD [2] is a simple realisation of such a discretised 
version of QCD and has been used for a long time in lattice simulations. 
However, it has been realised that this formulation has a number of severe 
problems: it shows large discretisation effects [3] that are linear in the lat- 
tice spacing a, it violates chiral symmetry strongly and develops unphysical 
small eigenvalues of the corresponding lattice Wilson-Dirac operator, even at 
rather large values of the quark mass. 

The problem of discretisation effects can be overcome following the Sy- 
manzik improvement program [4, 5] , introducing the well known clover term 
[6] to get full, on-shell 0(a) improved results, if the improvement is performed 
non-perturbatively [3]. In this approach then also the chiral properties are 
improved, although chiral symmetry breaking and 0(a 2 ) lattice spacing ef- 
fects are still left in the theory and have to be extrapolated away. Despite 
the fact that non-perturbatively improved Wilson fermions clearly diminish 
discretisation errors, they, unfortunately, show the same - if not worse - 
problem of the appearance of small eigenvalues of the lattice Dirac operator 
[7]- 

In order to solve the problem of small eigenvalues, it has been proposed 
in [8] to use the so-called twisted mass formulation of QCD. In this approach 
the mass term in the Dirac operator is chirally twisted [8] , see also [9] . When 
using such a chirally twisted lattice action in combination with the clover 
term and a non-perturbatively tuned value of the improvement coefficient, 
the theory is 0(a) improved and the corresponding lattice Dirac operator is 
safe against developing small eigenvalues. 

Recently, it has been realised [1] that all the above properties of non- 
perturbatively improved, twisted mass QCD can also be obtained when the 
clover term is completely omitted. If a special value of the twisting an- 
gle is chosen one theoretically obtains 0(a) improved results without adding 
improvement terms. At the same time, the lattice Dirac operator is still 
protected against the appearance of small eigenvalues by construction. This 
curious observation receives a special importance for simulations with dy- 
namical quarks: In [10] it was found that dynamical fermion simulations 
with non-perturbatively improved Wilson fermions show signs of first order 



2 



phase transitions which render simulations very difficult and induce large 
cut-off effects. 

Although by varying the form of the gauge actions [11] the first order 
phase transition seems to vanish, it is unclear, whether eventually such phe- 
nomena will reappear. The problem of the presence of the first order phase 
transition may be related to the fact that the clover term in the fermion 
action generates an adjoint gauge action. Therefore, if 0(a) improvement 
can be achieved without the clover term, as anticipated in [1], the problems 
connected to such phase transitions should be completely eliminated. As a 
consequence, the potential of twisted mass fermions in general may be very 
large. Since the small eigenvalues are regulated by the twisted mass pa- 
rameter, simulations at much lower quark masses than used today could be 
performed with promising, but hard to estimate, advantages to explore the 
chiral limit of lattice QCD. 

In this paper we provide a first test of the conjecture of 0(a) improve- 
ment of Wilson twisted mass QCD in quenched numerical simulations. For 
this purpose we performed a scaling test of the vector meson mass my and 
the pseudoscalar decay constant F PS at a fixed value of the physical pseu- 
doscalar mass mpg. We compare the results with those that have been ob- 
tained for standard Wilson fermions, see [12] and references therein, and 
non-perturbatively improved clover fermions [13]. 



2 Twisted mass QCD with Wilson quarks 

In twisted mass QCD (tmQCD) as formulated in [8] the twisted mass action 
in the continuum reads as follows: 

S F [ip, $ = J d 4 xi)(D^ il + m + ifiq^T 3 )^ , (2-1) 

where denotes the usual covariant derivative, m is the standard bare 
quark mass, r 3 is the third Pauli matrix acting in flavour space and [i q is the 
twisted mass parameter, also referred to as the twisted mass. 

An axial transformation, 

ip' = exp(zcJ7 5 r 3 /2)V>, ijj' = vjjexp(iu;>y 5 T 3 /2) , (2-2) 

with a real rotation angle uj leaves the form of the action invariant and merely 
changes the mass parameters into m' Q and fi' , 

m' = m cos(u;) + \i q sin(o;) 

/i' q = — m sin(cj) + ji q cos(u;) . (2-3) 
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The standard action (// = 0) is obtained by setting tancj = fi q /m . Note 
that t 3 is traceless and therefore the transformation (2-2) does not couple to 
the fermion determinant anomaly. 

In order to have an 0(a) improved twisted mass lattice action, it appears 
to be natural to discretise the Dirac operator adding appropriate improve- 
ment terms to the standard Wilson-Dirac operator. Indeed, it has been 
shown in [14, 15, 16, 17] that by adding the usual clover term, full on-shell 
0(a) improvement can be obtained. 

In [1] it has been realised later that using simply the standard lattice Dirac 
operator Dw one can obtain 0(a) improved physical observables without 
adding improvement terms. More precisely, it is possible to obtain 0(a) 
improved lattice results by employing only the standard massless Wilson- 
Dirac operator, 

Dw = ^{7m(V; + V„) - arV;V M } (2-4) 

under the condition that one averages over physical observables that are 
obtained from simulations at positive and negative values of the Wilson pa- 
rameter r. A less general, but similar suggestion was made by the authors 
of refs. [18, 19]. Instead of using positive and negative values of r, quark 
masses with different signs may be used: the bare quark mass mo can be 
written as 

m = m c {r) + m q , with m c (— r) = — m c (r) , (2-5) 

where m c (r) is the critical quark mass, and m q is the subtracted bare quark 
mass. Averaging physical observables obtained from simulations at positive 
and negative subtracted bare quark masses, again 0(a)-improvement is ob- 
tained. 

Let us shortly sketch the arguments leading to this surprising result. One 
first has to observe that with 

^5 : < j j, j (2-6) 

the following combined transformation 

Tlf = TZ 5 x[r^ -r] x [m q -> -m q ] (2-7) 

is a so called spurionic symmetry of the ordinary Wilson action. Another 
symmetry of the Wilson (and Wilson tmQCD) action is TZ$ x T> d . In the 
continuum, the transformation T>d has the effect of changing the sign of all the 
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space-time coordinates and multiplies each local term £j in the Lagrangian 
density by the factor (— \) d ^ Ci \ where d[Ci] is the naive dimension of £j. The 
lattice version of this transformation is more involved and we refer to ref. [1] 
for details. 

Taking now the parity properties of multiplicatively renormalisable op- 
erators under TZ^ and 71% x T>d into account, one can show - using the 
Symanzik expansion - that one gets 0(a) improvement when averaging over 
two simulations with positive and negative Wilson coefficient r (Wilson av- 
erage (W^4)). In addition, from the spurionic symmetry TZ S ^ of the Wilson 
(and Wilson tm) action one can obtain 0(a) improved physical observables 
when averaging, at a fixed value of r, over two simulations with positive and 
negative m q , as defined in eq. (2-5) (mass average (MA)), taking into account 
this time the ^-parity of the operators. Studying the chiral properties of 
the scalar condensate with Wilson fermions, a similar suggestion was made 
by the authors in ref. [20]. In the special case of choosing u> = ±7r/2, such 
an averaging procedure is done automatically. A change of the sign of r 
is equivalent to uj — > uj + n. Hence, for uj = ±7r/2 all the quantities that 
are even under to — > —uj are automatically improved without any averaging 
procedure. 

It is the aim of this paper to check this conjecture in practical simulations 
by performing a scaling test for the vector meson mass my and the pseu- 
doscalar decay constant Fpg at uj = tt/2. The main goal is to test whether 
the results for m v and F PS are consistent with the anticipated leading 0(a 2 ) 
behaviour and that the linear a dependence is indeed cancelled. In addition, 
it is an interesting and important question, what the size of the remaining 
lattice spacing effects arising in 0(a 2 ) will be. 

Let us list a few properties of the composite fields in the twisted mass 
formulation before going to the numerical results. Due to the transformation 
rule (2-2) one also has to transform the composite fields defined in the usual 

way, 

S°(x) = <Kx)V(x), r(i) = ^W7 5 ^(4 

= ^(x) 7M 75y^(x), V«(x)=?P(xh^(x). (2-8) 

As an example we give here the relations for the axial and vector currents 
in the "physical basis" (primed quantities) and "the twisted basis" (unprimed 
quantities), 

_ jcos(uj)A« + e 3 ^ sm{u>)Vg (a = 1,2), 
"~U (« = 3), ( } 
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V" 



la 



cos{u)V* + e 3af3 sin(o;)Aj (a = 1, 2), 
V; 3 (a = 3). 



(2-10) 



Note that in eq. (2-9) for a = 1,2 and u = tt/2 the role of the axial and 
vector currents are just interchanged. Of particular interest is the PCVC 
relation, which takes the following form in the twisted basis: 



where d* is the usual backward derivative. Through a vector variation of the 
action one obtains the point-split vector current as defined in [8, 16]. This 
current is protected against renormalisation and using the point-split vector 
current, the PCVC relation is an exact lattice identity. This implies that 
Zp = Z~ x , where is the renormalisation constant for the twisted mass 
fi q . This will become important in the extraction of the pseudoscalar decay 
constant Fps as described below. 

3 Numerical Tests 

In this section, we describe our numerical results for testing the scaling 
behaviour of Wilson tmQCD in the quenched approximation. We started 
our investigation by performing a comparative benchmark study of different 
solvers for obtaining the quark propagator. We found the CGS algorithm 
[21] to be superior to the BiCGstab and the CG algorithms. We therefore 
used the CGS algorithm throughout this work. Gauge field configurations 
were generated by standard heat-bath and over-relaxation techniques. 

3.1 Mass average 

In order to test the predictions of ref. [1], we started with the mass average 
procedure. To this end, we selected a value of j3 — 5.85, set the Wilson 
parameter r = 1 and performed simulations on 12 3 x 24 lattices at positive 
and negative values of m q = 1/2 (1/k— 1/k c ). While for m q = +0.02725, 
the propagator computations went smoothly, for m q = —0.02725 the com- 
putation of the quark propagator was exceedingly expensive. The reason for 
this behaviour can be traced back to the spectrum of the Wilson tmQCD 
operator as can be seen in fig. 1. 

Comparing fig. 1(a) with fig. 1(b), in the case of negative m q one has to 
deal with extremely small eigenvalues of the operator (LV + m o)^(£V + mo)- 
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-2/2 q e 3aP P f3 



(2-11) 
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Figure 1: Monte Carlo time evolution of the eleven smallest eigenvalues A of (Dw + 
mo) (Dw + m o), normalised by the largest eigenvalue, at fj. q — and m q = ±0.02725 on 
a 12 3 x 24 lattice (k c = 0.161662(17), = 5.85). 
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Clearly, these very small low-lying eigenvalues lead to a poor convergence of 
the solver employed and hence to very costly simulations. Projecting out 
these small modes does also not help in this respect since the computation 
of the eigenvalues is again costly 

We therefore proceeded to the "self- averaging" case of choosing uj = tt/2, 
for which it has been shown in [1] that one gets 0(a) improvement even 
without the need of any averaging for all quantities that are even under 
uj — > —uj. In our practical implementation we have used the twisted basis. 
Hence, the choice uj = 71/ 2 corresponds to set m q = and « g ^ 0. In this 
situation, the corresponding Wilson-Dirac operator is protected against small 
eigenvalues and we do not expect difficulties with the simulations. 

3.2 Scaling of the vector meson mass 

In order to verify the prediction of ref . [1] we computed the vector meson mass 
my and the pseudoscalar decay constant Fp$ for the following values of j3: 
5.85, 6.0, 6.1, 6.2. We used periodic boundary conditions. The corresponding 
lattice volumes were 14 3 x 28, 16 3 x 32, 20 3 x 40 and 24 3 x 48, respectively. 

For our simulation in the twisted basis at m = m c we had to determine 
the critical hopping parameter k c for each value of (3. At all the (3 values of 
our simulations we made our own determination of the value of n c from the 
intercept in n at zero pion mass. The values of k c are given in table 1. Note 
that these critical values of k have an intrinsic uncertainty of 0(a). This 
is, however, sufficient for obtaining fully O(a) improved results in Wilson 
tmQCD [1]. 



p 


L 


T 






N 

1 ' meas 


5.85 


14 


28 


0.161662(17) 


0.0376 


400 


6.0 


16 


32 


0.156911(35) 


0.03 


388 


6.1 


20 


40 


0.154876(10) 


0.025854 


299 


6.2 


24 


48 


0.153199(16) 


0.021649 


215 



Table 1: Parameters of the simulations. Note that the values for k c are obtained from a 
different set of measurements. 

In order to fix the physical situation in our scaling test, we kept r mpg 
fixed for all values of f3. For this purpose we determined the value of [i q to 
fix romps = 1.79. The corresponding values of [i q for each value of f3 and all 
our simulation parameters are given in table 1. 

We computed the standard 2-point correlation functions at zero momen- 
tum for the pseudoscalar and axial operators (which in the twisted basis at 
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uj = 7r/2 gives the correct operator to extract the vector meson mass), 

fp(t) = X>°(^(0)> (3-1) 

X 

/IW^EEWMP)) (3-2) 

i=l x 

with P a (x) and A a (x) given in eqs. (2-8) and x = (x, t). In order to obtain a 
non- vanishing result, the flavour index has to be the same in these correlation 
functions and we will choose a = 1 in the following. The pion mass could 
be extracted easily from the exponential decay of the correlation function 
fp{t). For the vector meson mass, we performed two mass fits for the ground 
state mass and the first excited state. We checked the stability of the fit by 
changing the value of t m i n where the fit started. As a cross check, we also 
determined the effective ground state mass and found consistent results. All 
errors were computed by a jackknife analysis. The numerical results at our 
simulation points are collected in table 2. 



p 


r /a 


ampg 


amy 


aFp S 


5.85 


4.067 


0.4340(16) 


0.656(11) 


0.1147(11) 


6.0 


5.368 


0.3329(21) 


0.488(11) 


0.0859(9) 


6.1 


6.324 


0.2871(17) 


0.427(9) 


0.0717(8) 


6.2 


7.360 


0.2438(16) 


0.363(10) 


0.0640(10) 



Table 2: Results for the vector meson mass and the pseudoscalar decay constant. 

In fig. 2 we show our results for the vector meson mass as a function 
of a 2 represented by the open circles. In addition, we also show results 
from non-perturbatively O(a) improved Wilson fermions (filled circles) [13]. 
Finally, we added results for standard Wilson fermion simulations, see [12] 
and references therein, as they were available in the literature. We remark 
that the published data were not always at exactly the same value of romps 
that we used for our Wilson tmQCD simulations. In such cases we performed 
an interpolation to the desired value of romps- The error from this (small) 
interpolation is negligible for the results presented here. 

We performed a simple extrapolation of the Wilson tmQCD and the O(a) 
improved data of the form r Q my = r my nt + b- (a/r ) 2 , with r my nt the con- 
tinuum value of the vector meson mass and r — 0.5 fm. For the pure Wilson 
results we replaced the quadratic term with a term proportional to (a/r ) in 
the extrapolation. Let us remark that our data for my for Wilson tmQCD 
has about a factor of four less statistics than the data from 0(a)-improved 
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2.0 ' 

0.00 0.02 0.04 0.06 0.08 

(a/r ) 2 

Figure 2: Scaling behaviour of the vector meson mass as a function of the lattice spacing 
squared for fixed pion mass, rgmps = 1.79. Open circles denote Wilson tmQCD, while 
filled circles are from non-perturbatively improved Wilson fermions [13]. The Wilson data 
without improvement (open squares) are collected from several sources in the literature. 
Note that the filled circles are slightly displaced for better visibility. 

Wilson fermions, which is reflected in the larger error bars. Nevertheless, it 
is evident that the Wilson tmQCD results show a very similar scaling be- 
haviour as the 0(a)-improved Wilson fermions. This becomes even clearer 
when we compare with the unimproved pure Wilson data for my that we 
show as open squares in fig. 2. Here large lattice artefacts are seen and the 
scaling behaviour is much worse than with Wilson tmQCD or 0(a)-improved 
Wilson fermions. We also note that the data for Wilson tmQCD are rather 
flat as a function of a 2 indicating that also higher order lattice spacing effects 
are suppressed. Clearly, it would be desirable to test these promising results 
in more precise simulations using a much higher statistics. 
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3.3 Scaling test for Fps 



In order to extract the pion decay constant Fps as another quantity to test 
scaling, we start with the standard definition of Fp$ (again fixing the flavour 
index a = 1), 

(0\A l \PS) =m PS Fps . (3-3) 

In the twisted basis, at uj = n/2, the axial current is related to the vector 
current by the transformation eq. (2-2), see eqs. (2-9,2-10), and so we can 
write 

^<0|V^ 2 |PS> = F PS m 2 s . (3-4) 

Using the vector Ward identity in eq. (2-11), we can finally relate the diver- 
gence of the vector current to the pseudoscalar density and obtain 

Fpsmps = <9 M (0|V; 2 |PS> = 2 f i q (0\P 1 \PS) . (3-5) 

For asymptotic Euclidean times, the pseudoscalar correlation function fp 
assumes the form 

/i ( t ) = l(0| f 1|PS)|2 . (e" m ^ + e -mps(T-m , a « t « T . (3-6) 
2m PS 

Thus, by fitting the pseudoscalar correlation function for large time separa- 
tions, we can obtain mpg and the amplitude | (0|P 1 |PS) | 2 / m PS from which 
we then compute the desired matrix element |(0|P 1 |PS)|. Hence, we have 
all necessary ingredients to determine Fp$ from eq. (3-5), without the need 
of any renormalisation factor, since Zp = Z" 1 , as we have mentioned in the 
previous section. The error estimate of the so computed value of F PS is per- 
formed with a jackknife procedure. We note that the discussion of how to 
get Fps with tmQCD resembles very closely the strategy one would follow in 
the continuum or with overlap fermions. 

In fig. 3 we show our results for ro-Fps at a fixed value of romps = 1.79. We 
also add results from 0(a)-improved Wilson fermions [13] in the figure. We 
performed an extrapolation of the Wilson tmQCD and the 0(a)-improved 
data of the form ro-Fps = Hj-Fpg* + b' ■ (a/ro) 2 , with roF'pg the contin- 
uum value of the pion decay constant in units of r ~ 0.5 fm. While the 
0(a)-improved Wilson results do show cut-off effects, presumably due to the 
particular value of the improvement coefficient ca employed in these calcula- 
tions, the results for Wilson tmQCD are essentially flat, indicating that also 
the higher order lattice spacing effects are substantially smaller. 
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Figure 3: Scaling behaviour of the pion decay constant Fps as a function of the lattice 
spacing for fixed pion mass, romps = 1.79. Open circles denote Wilson tmQCD, while 
filled circles are from non-perturbativcly improved Wilson fcrmions [13], for which only 
the three data points with smallest values of (a/ro) 2 are included in the fit. Note that 
the points corresponding to the continuum extrapolation are slightly displaced for 

better visibility. 
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4 Conclusion and Outlook 



In this paper we tested the striking idea of using the standard Wilson-Dirac 
operator with a twisted mass and twist angle uj = it/ 2 to obtain full 0(a) 
improvement for physical quantities that are even under to — > —uj. The 
results of our - quenched - study for the vector meson mass my and the 
pion decay constant Fps are very encouraging. It seems that in this setup 
the lattice spacing effects are substantially reduced with respect to standard 
Wilson fermions and consistent with vanishing 0(a) discretisation errors. At 
the same time it seems that also lattice spacing effects that come in higher 
orders in a are small, a result that could not be anticipated before. 

Thus it seems that Wilson twisted mass QCD has the potential to solve 
many problems at once: 

• It automatically reduces lattice spacing effects, consistent with 0(a) 
improvement. 

• It protects against very small low-lying eigenvalues of the lattice Dirac 
operator. 

• It can avoid unwanted phase transitions as appeared in dynamical simu- 
lations with clover improved Wilson fermions and Wilson gauge action. 

In the light of this, dynamical simulations with twisted mass QCD ap- 
pear to be very promising to approach the physical point at realistic values 
of quark masses. In some way, the twisted mass plays a similar role as an 
infrared cut-off as the quark mass in the staggered fermion or overlap ap- 
proach does. The very good performance of staggered dynamical simulations 
gives hope that also in twisted mass QCD the simulations can be acceler- 
ated as compared to standard, improved or unimproved, Wilson fermions. 
Naturally, with Wilson twisted mass QCD the problem of unphysical taste 
degrees of freedom is completely avoided, which, to us, is a clear advantage 
of the tmQCD idea. 
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